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<N : Abstract 

> : 

^sO . Gauge theories broken by a single Higgs field are known to have first-order 

(N ■ 

^ ! phase transitions in temperature \i \/ <^ 1, where g is the gauge cou- 

^ ! 

\^ , pling and A the Higgs self-coupling. If the theory is extended from one to 
0\ ; 

i N Higgs doublets, with U(A^) fiavor symmetry, the transition is known to be 

O-t'. 

^ ', second order for X/g'^ ^ 1 in the N ^ oo limit. We show that one can in 

(D ! 

r~| , principal compute the tricritical value of \/g , separating first from second- 

^ ■ 

I order transitions, to any order in In particular, scalar fluctuations at 

^ , the transition damp away the usual problems with the infrared behavior of 

high-temperature non-Abelian gauge theories. We explicitly compute the tri- 
critical value of X/g'^ for U(l) and SU(2) gauge theory to next-to-leading order 
in 
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I. INTRODUCTION 



The possibility of explaining the observed baryon number of the universe by physics 
occurring at the electroweak phase transition has, in recent years, renewed interest in un- 
derstanding how to assess the existence, order, and strength of phase transitions in gauge 
theories, such as electroweak theory, where gauge bosons get mass by the Higgs mechanism. 
The phase structure of the electroweak sector of the minimal standard model, with a single 
Higgs doublet, makes a good starting point for exploring such transitions. It has long been 
appreciated []T| that the phase transition is first-order in the limit that the zero-temperature 
Higgs boson mass is small compared to the W boson mass, i.e. when A <^ g"^, where g is the 
electroweak gauge coupling and A is the Higgs self-coupling. In this limit, a perturbative 
analysis of the phase transition is adequate to establish its order and compute its physical 
properties. (Here and throughout, we assume A and g^ are both small.) What has been 
more difficult is to study the transition when X ^ g^. In this limit, perturbation theory 
breaks down due to large infrared fluctuations characteristic of critical or near-critical be- 
havior. One response is to turn to numerical simulations of the transition. It is interesting, 
however, to see what can be said about the transition analytically if one modifies the theory 
to make it more tractable. For example, if the three spatial dimensions are replaced by 
4— e dimensions, where e <^ 1, then it is known that the transition remains first-order for 
any finite X/g"^ 10)131 • If the Higgs sector is generalized to contain Higgs doublets with 
U(A^) symmetry, then in the N ^ oo limit the transition is first-order for X/g"^ <^ 1/N 
and second-order for X/g"^ ^ |3|- Recent numerical simulations for A^=l, in contrast, 
suggest that the first-order transitions end at a critical value of X/g^ above which there is 
no phase transition whatsoever 



^ We should emphasize that one of the pieces of evidence presented in ref. — the large volume dependence 
of the susceptibility — has a loophole which is realized in a simple and relevant example. Ref. |^ shows 
evidence that this susceptibility approaches a constant in the large volume limit, for Higgs masses above 
roughly 80 GeV. In a second-order transition, the large volume behavior of this susceptibility should be 
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The goal of the present work is to extend understanding of the large limit beyond 
leading order in studying in particular the critical value of \/ demarking the end of 
first-order transitions. We emphasize that large here refers to the number of scalar fields 
and not to the replacement of the gauge group by SU(A^). 

Usually, studying critical behavior of weakly coupled field theories is more difficult than 
studying those theories far from the transition, because long-distance fluctuations appear 
at the transition whose physics is non-perturbative. Small pure scalar theories, for ex- 
ample, are easy to study far from the transition but difficult near the transition for this 
reason. Amusingly, large A^(scalar) non-Abelian gauge theories are exactly the opposite: At 
temperatures far above the transition, electric forces are Debye screened in the hot plasma 
but magnetic forces are not, and the non-Abelian nature of the forces gives rise to mag- 
netic confinement at large distances, which cannot be treated perturbatively. But, as we 
shall discuss, at the phase transition long-distance scalar fluctuations screen the magnetic 
forces sufficiently to prevent magnetic confinement. The long-distance scalar fluctuations 
themselves are treatable in a 1/A^ expansion just as in large A^ pure scalar theories. 

In the remainder of this introduction, we briefly discuss how large one might suspect A^ 
has to be for the large A^ expansion to be useful. Then we discuss whether a moderately 
large A^ Higgs sector is phenomenologically viable. In section 2, we will flx notation and 
briefly review that the problem of flnite temperature phase transitions in 3-1-1 dimensions 
is equivalent to the study of fleld theories in 3 Euclidean dimensions. Then we discuss the 
power counting of the loop expansion and why magnetic conflnement is not a problem at 



an analytic function of V^^ plus a non-analytic scaling piece whose leading term is V"'/^^ ^ where a and 
v are respectively the specific heat and correlation length exponents. However, there are systems with 
second-order transitions where a/v is negative^ and so the measured susceptibility would indeed approach a 
constant rather than diverge. A relevant example is the pure scalar sector of electroweak theory itself, the 
0(4) model, where a/iv = —0.33(4) (see table |). A more convincing demonstration of the absence of a 
transition is the plot in ref. Q of the inverse correlation length vs. temperature, which shows no suggestion 
of a divergence in the correlation length. 
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the transition. We will also see that the calculation of the tricritical value of X/g"^ order 
by order in is conceptually more straightforward than the calculation of many other 
quantities. In section 3, we carry out this computation to next-to-leading order for the U(l) 
gauge theory. Section 4 is devoted to clearing up some minor subtleties of regularization of 
diagrams. Finally, we carry out the next-to-leading order computation for SU(2) theory in 
section 5. same for SU(2) theory in section 5. 

A. How large in large A^? 

If gauge interactions are ignored, the scalar sector of the minimal standard model with 
a single Higgs doublet is equivalent to an 0(4) theory of four real scalar fields and has a 
second-order transition. Table | demonstrates large results applied to this case ((7=0). 
Amusingly, next-to- leading order in actually gets in the right ballpark for various critical 
exponents. Sadly, this happy circumstance that 1 is just on the verge of being a large number 
of doublets will not survive the inclusion of gauge interactions. 







Large 




actual 




LO 


NLO 


NNLO 


series 


monte carlo 


7 = 


2 


1.392 


1.188 


1.44(4) 


1.477(18) 


V = 


1 


0.730 


0.612 


0.73(2) 


0.7479(90) 


/3 = 


0.5 


0.399 


0.325 


0.38(1) 


0.3836(46) 


6 = 


5 


4.594 


4.695 


4.82(5) 


4.851(22) 


V = 


0.0675 


0.0554 


0.0260 


0.03(1) 


0.0254(38) 


a = 


-1 


-0.189 


0.163 


-0.19(6) 


-0.244(27) 



TABLE I. Large N expansion results [|| for critical exponents in the pure scalar case (5=0) for one 
Higgs doublet: the 0(4) model. Results are given for leading (LO), next to leading (NLO), and next to 
next to leading (NNLO) order in 1/N. The actual values, as predicted from series analysis ( |^ as cited in 
[0) or measured by monte carlo are also shown. There are various scaling relationships between these 
quantities, so that only two of the above exponents are independent. 
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At another extreme, one can analyze the phase structure for arbitrary A^ in 4 — e spatial 
dimensions [^,|^. One finds that the qualitative picture given by the large A^ limit — that 
there is a tricritical value of X/g'^ above which the transition is second order — is correct 
when 

N > 182.95 - 320. 50e + 0(6^) , for U(l) with N charged scalars , (1.1) 
A^ > 359 - 495.4e + 0(6^) , for SU(2) with A^ scalar doublets . (1.2) 

So, near four spatial dimensions, large A^ means A^ ^ 183 or A^ ^ 359, respectively! 

A result of our calculation of the tricritical point will be that quantitative success of 
large A^ in three dimensions appears to require A^ ^ 4 charged scalars for U(l) theory and 
A^ > 20 doublets for SU(2) theory. 

B. Is moderately large N phenomenologically viable? 

Generalizing the one Higgs model to a U{N) Higgs model is motivated solely by the 
desire to find a theory whose phase transition is analytically tractable. Nonetheless, it's 
interesting to briefly consider (just for fun!) whether such a model might be consistent with 
real world phenomenology. 

A U{N) flavor-symmetric Higgs sector is a phenomenological disaster because elec- 
troweak symmetry breaking will break the global U{N) and produce massless Goldstone 
bosons. The essential difference between Higgs models susceptible to large A^(scalar) analy- 
sis and generic multiple Higgs models is the necessity of a large A^ global symmetry. There 
is no reason, however, that this symmetry need be continuous. Though we do not study it 
in this paper, one could instead have a Higgs sector with a discrete A^-flavor permutation 
symmetry:^ 



^ The purely scalar sector is a simple generalization of the cubic anisotropy model. See, for example, 
ref. 0. 
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i 

^9uQ^{^X + • ■ ■ + $7v)ml + ^dOR^2($l + ■ ■ ■ + $7V)*4 • (1.3) 

As an added bonus, the permutation symmetry prevents the tree-level flavor- changing neu- 
tral currents that plague generic multiple scalar models. 

Because the discrete flavor symmetry is spontaneously broken, the model suggested above 
will produce cosmological domain walls which overdose the universe. This problem could 
by solved by the introduction of a very small symmetry breaking term (which is natural in 
the sense of 't Hooft [jrT|) that would cause the domain walls to coalesce after they were 
formed. 

Some sort of > 1 models therefore seem acceptable phenomenologically. It is worth 
noting that there is an important qualitative difference between the iV = 1 and > 1 cases. 
For A^ > 1 there is, by construction, a global flavor symmetry that is spontaneously broken 
at the electroweak scale. Such models therefore always have some sort of phase transition 
(ignoring the tiny symmetry breaking term discussed above). The A^ = 1 model, in contrast, 
need not have a transition because, technically, local symmetries are never spontaneously 
broken due to Elitzur's theorem 

We shall not analyze in any detail just how large A^ could be in a realistic theory except to 
make one, trivial observation: there is a simple constraint from triviality. Non-perturbative 
continuum scalar theories are not well deflned, and the effective strength of scalar interactions 
is 0{NX) instead of 0(A). The largest Higgs mass for which the theory can be sensible as 
an effective theory therefore decreases roughly as from the A^=l limit of 0(1 TeV). 

A^ ^ 100 is clearly out of the picture. On a related note, A^ > 20 would destroy the 
asymptotic freedom of the SU(2) electroweak interactions in the standard model. 

One can imagine that A^ might be big enough for the large A^ approximation to be 
reasonable, but not too big to run into phenomenological problems. There's no good reason, 
of course, why nature would choose to be so peculiar. 
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II. LARGE N COUNTING 



The problem of studying a second-order (or very weakly first-order) phase transition in 
weakly coupled quantum field theory, as one varies the temperature, can be reduced to the 
problem of studying the phase transition in three-dimensional Euclidean field theory, as one 
varies a mass. In our case, the three-dimensional theory is of the form 

S = J d^x{\D$\^ + + m^\$\^ + iA|$|^ + ^r/|$r} , (2.1) 

and the transition occurs as is varied through zero. Let's take a moment to briefly review 
this correspondence. 

Critical or near-critical behavior of phase transitions is governed by the physics of long 
distances, which in our case is classical. It is classical because the Bose density per mode, 
l/(e~^^ — 1), becomes large for small energies E. So one may study the long-distance, 
equilibrium properties of such transitions by studying the classical statistical mechanics of 
field theory in three spatial dimensions. More formally, one may start with the Euclidean 
formulation of finite-temperature quantum field theory and then integrate out the physics 



of the small, periodic Euclidean time direction. (See refs. for a review.) Taking the 
additional step of integrating out the Debye-screened Aq field, one obtains an effective three- 
dimensional theory of the form ( p2.1| ) plus irrelevant interactions. The parameters of the 
effective three-dimensional theory can be perturbatively related to those of the original 3-1-1 
dimensional theory if the couplings are small. One finds, up to higher-order corrections, 
that 

g' = gl{T)T, X = X,iT)T, r] = 0{glX'), (2.2) 

where gi{T) and \a{T) are the dimensionless couplings of the 3-1-1 dimensional theory at 
a renormalization scale of order the temperature T. One also finds that increasing the 
temperature through the transition corresponds to varying the scalar mass m in the effective 
theory from < to m? > 0. So the problem of understanding the phase transition of the 
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original theory in temperature is equivalent to understanding the phase transition of a three- 
dimensional theory in m^. In the limit that the original couplings gf and A4 are considered 
arbitrarily small, the problem of finding the tricritical value of Xa/qI in the original theory 
is the same as finding the tricritical value of X/g'^ in the three-dimensional theory. 

So focus on the three-dimensional theory ( ^2.11 ), and note that A and have dimensions of 
mass. Now consider the naive perturbative expansion for some physical observable associated 
with a small momentum scale p and suppose we are in the symmetric phase > and that 
p <^m. Then the scalars will decouple, and we must focus on the non-Abelian interactions 
of the magnetic gauge fields. (As mentioned earlier, the electric ones are Debye screened.) 
By dimensional analysis, the loop expansion parameter for gauge interactions is then g'^/p, 
and perturbation theory will fail once we try to explore momentum scales p g"^. This is 
the source of the infrared problem for non-Abelian gauge theories at high temperature. 




FIG. 1. Leading contribution to the gauge boson self-energy. 

Now consider the case right at a second-order transition, where the scalar mass is zero, 
and consider the effect of the self-energy diagrams of fig. |l] on the gauge boson propagator. 
By dimensional analysis, this self-energy is 



T{{p) = aNg^p, 

where a is a numerical constant, and the gauge propagator becomes 

1 1 



(2.3) 



G{p) 



as p 



0. 



(2.4) 



p"^ + aNg^p aNg^p 
This is less divergent in the infrared than the perturbative propagator l/p"^. For distances 
r ^ l/Ng"^, the scalar degrees of freedom have screened the gauge propagator from 1/r 
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behavior to behavior. The problematical interactions above the phase transition were 
non-Abelian gauge interactions. Now, with the propagator ( p.4| ), such interactions are under 
perturbative control for large A^. Loops of gauge bosons will in general be infrared conver- 
gent, and the scale of the loop momenta will be 0{Ng'^) if the external momentum is small. 
The cost of adding a new pair of non-Abelian interactions to a graph, such as depicted in 
fig. @, is then, by dimensional analysis, generically 



The double lines represent the resummed gauge propagator ( p.4|) . 



(2.5) 





FIG. 2. Adding a pair of non-Abelian interactions to a graph. 

We have now discussed the cost of adding purely gauge loops to a diagram. Before 
proceeding to the case of generic scalar loops in a diagram, it will be useful to first review 
the leading-order calculation of ref. ||^ for the tricritical value of X/g"^. 

A. Tricritical point at leading order 

For a second-order phase transition, the transition occurs when the effective mass of the 
scalar field vanishes. A tricritical point occurs when the effective, low-momentum, quartic 
coupling Aeflf of the scalar vanishes as well. The classic mean-field example is the potential 



(2.6) 



10 



which (ignoring corrections due to fluctuations) has a second-order transition in if A > 0, 
a first-order transition in if A < 0, and a tricritical point at A = 0. The actual low- 
momentum effective potential for gauge-Higgs theories was computed at leading order in 
in ref. [0]. Here, we just need the effective value of the four-point interaction. We 
begin by assuming \/ g"^ is 0(A^^^), which we shall see a posteriori is the correct place to 
look for the tricritical point. For simplicity, we will also ignore the bare six-point coupling 
rj by setting it to zero.0 The theory is then super-renormalizeable. A and do not require 
renormalization and will henceforth refer to their bare, short- distance values. 

It is convenient to henceforth think of g"^ as 0{N~^) and so A as 0{N^'^). This is just 
a convention because g^ is dimensionful, but it is a convenient convention because it makes 
the internal momenta of gauge bosons in the graphs discussed above 0{N^). N counting of 
those graphs then reduces to counting scalar loops and explicit coupling constants. 



X + 




FIG. 3. Leading-order contributions to the scalar four-point interaction. 



The leading-order graphs for the four-point interaction are shown in fig. ^, where the 
double lines again represent the resummed gauge propagator (P^). For the U(l) case, these 
graphs give 

d?p 1 Qg'^ 



Ks = \-I2g^ j 



A 



(27r)3 (p2 + aNg'^pY n'^aN 
where computation of the self-energy diagrams fig. |l| yields 



(2.7) 



^ This is sensible, based on (|]^), if one formally considers gl to be arbitrarily small. If one instead considers 
the natural choice that is 0{N^^), then rj will be 0{N~^). At the order under consideration, the effects 
of 77 on the following derivation can be absorbed into the definition of A. 
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(2.8) 



Setting Aeflf to zero, the tricritical point is at 



A 96 



+ 0{N-^). 



(2.9) 



g2 T^2fq 



The case of SU(2) with doublet Higgs bosons differs just by the number of gauge bosons 
and the normahzation of the couphng of the scalars to the gauge bosons. With conventional 
normalization of 



Consider for a moment pure scalar theory. At zero external momentum, the naive loop 
expansion parameter is NX/m, where the 1/m follows from dimensional analysis. There 
is therefore an infrared problem when one goes to the transition, m 0. The standard 
application of large N techniques to this theory eradicates this problem by a large N resum- 
mation of the quartic interaction, which curbs the infrared behavior of the diagrammatic 
expansion. 

In gauge-Higgs theories, the same problem potentially arises and will require a similar 
resummation of the scalar interactions. However, for the special case of computing the 
tricritical value of X/g"^, this resummation is unnecessary. Because the problem with the 
naive expansion in scalar loops was an infrared problem, the loop expansion parameter 
NX/m due to the infrared behavior (p ^ Ng"^) of loops should be replaced by NX^^/m. 
But Acff at the tricritical point is zero by definition. As we have reviewed, that zero occurs 
because of a cancelation of interactions, such as those shown at leading order in fig. ^. This 



* This is most easily achieved by the standard technique of replacing the quartic interaction A(/)^ by 
+ \f\(^x where x is an auxiliary field, integrating out 0, and then studying the resulting theory of x- 



A 36 



+ 0(Ar-2). 



(2.10) 



g^ Ti^N 



B. Scalar loops in the infrared 
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cancelation breaks down for loop momenta p ^ Ng'^, so the real cost of adding a scalar loop 
at the tricritical point will be determined by the scale Ng"^: 



NX A 1 



(2.11) 




(a) 




(b) 




(c) 











(i) 







(k) 



FIG. 4. Next-to-leading order contributions to the scalar four-point interaction: abelian graphs. We 
have neglected graphs that vanish in the Abelian case due to Furry's theorem, i.e. charge conjugation. See 
fig. ^for interpretation of (f-k). 
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(1) 



(m) 



(n) 




(o) (p) 

FIG. 5. Next-to-leading order contributions to the scalar four-point interaction: SU(2) graphs in Landau 
gauge. We have neglected graphs which vanish by charge conjugation. Note that the subgraph of fig. |^ 
vanishes for SU(2) but not other groups. The dashed lines represent ghosts. See fig. ^ for interpretation of 
(l-n,p). 




FIG. 6. A subgraph which vanishes for U(l) and SU(2). 
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FIG. 7. Meaning of abbreviated graphs in figs. ^ and |[ 




FIG. 8. Example of diagrams which vanish in Landau gauge for zero external momentum. 

The upshot is that, to determine the tricritical point, we can set m=0 and then use 
naive large power-counting of diagrams, treating ~ 0{N^^) and A ~ 0(iV^^). The 
relevant diagrams for computing Agg at next-to- leading order in are shown in figs. ^ 
and ^ for Landau gauge. In this gauge, diagrams of the form of fig. § vanish at zero 
external momentum. The graphs (f-n,p) denote the contributions to Acg in the short-hand 
style explained by fig. |^. At each order in the calculation of Aeg, and the simultaneous 
determination of the tricritical point by setting Aefr = 0, one will find that scalar infrared 
divergences always cancel at the order one is calculating. We now turn to a calculation of 
the U(l) case, corresponding to the diagrams fig. ^, where this will be made explicit. 
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III. U(l) TRICRITICAL POINT: NEXT-TO-LEADING ORDER 



For simplicity, we shall present the calculation in Feynman gauge. The total contribution 
of diagrams of the form of fig. ^ still vanishes in the U(l) case because of the U(l) Ward 
identity of fig. (Landau gauge results, which will be useful for the SU(2) case, are given 
in Appendix |B 2|.) 




FIG. 9. A U(l) Ward identity. 

We shall proceed by doing the scalar loop integrals of all the diagrams of fig. ^, summing 
up the diagrams, and then integrating over the gauge boson momenta. The order of the 
last two steps is important because the gauge momentum integrals are infrared divergent 
for the individual diagrams but not for the sum, and it will be convenient not to have to 
introduce a consistent infrared regulator. (Our disregard of regularization will sometimes 
be a bit cavalier in this section, and we delay discussion of potential subtleties to section 4.) 

As a warmup, consider fig. This diagram gives a contribution to Acs of 

a« = -72Ng' f T,, f \ = -Ng' f -^^^^ , (3-1) 

Jpq J«/^(l + p + q)^ Jpq |p + q| 

where we have introduced the notation 



p 



^pq = fpfq ~\~ fpfq + fpfq ; (3-3) 



and fp is the large N resummed gauge propagator 



^ We have checked this exphcitly. We have also exphcitly checked our final U(l) results are the same in 
any covariant gauge. 
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1 

+ aNg'^p 



fp ^ ZTTT^ ■ (3-4) 



A useful table of various / integrals is given in Appendix 0. 

As a slightly more complicated example, consider fig. ^, which gives 

aS^4B./^^./ g;,4.PV^ . (3.5, 

To simplify the / integration, one may use the standard technique of rewriting numerators 
in terms of denominators and things that don't involve /: 

(21-p)-(21 + q) = (l-pf + (l + qf + 2/2-(p2 + p-q + g2), (3.6) 

giving 

Jva J I 



^ 1 ^ 2 (p2 + p.q + g2) 



Jpq 



6 6 12 6(p2 + p.q + g2^ 



+ - + 



■ ■ , , , , (3-7) 

_q p |p + q| M|P + q| 

Appendix |^ explains an amusingly simple method for evaluating the last I integral by using 

a simple change of variables. 

Fig. can be done similarly, and the result is given in Appendix 0. 

Figs, ^(a-h) are slightly more subtle because the scalar integration is infrared divergent, 
both diagram by diagram and collectively. For instance fig. ^a gives 

(3.8) 

However, as we discussed earlier, A should end up replaced by Aefr = in the infrared if 
we sum up diagrams, as was shown at leading-order in fig. 0. The cancelation in fig. ^ will 
correspond to a cancelation, at this order in 1/A^, in the diagrams of fig. ^(a-c,f-h) if we 
consider the pieces of (f-h) represented by the second term on the right-hand side of fig. 0a.0 
That is, if we set 



^ Keep in mind that the last two terms of fig. || can be ignored on the context of the NLO diagrams (a-h) 
because of the Ward identity of fig. 
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A 



(3.9) 



to make Aefr zero at leading order, then we will find that the infrared divergences just 
discussed will cancel each other at the order in at which we are computing. It will be 
convenient to write this condition at a more primitive level, related directly to the diagrams 
of fig. Bl as 



X = 12C,g' f /|, 
Jp 

where the value of A is now parametrized by C\ and 



Cx = l + OiN 



at the tricritical point. 



(3.10) 



(3.11) 




+ 



+ 




+ 



— X— 



= 



FIG. 10. Diagrams contributing to the effective scalar mass rrics. 



FIG. 11. IR divergent loops caused by treating a non-zero mass perturbatively. The circles represent 
mass insertions on a massless propagator. 

There still remain infrared divergences, however, in the pieces of figs, ^(d-h) correspond- 
ing to the first term on the right-hand side of fig. |^. The divergence does not cancel in 
the sum of these graphs because, though we have taken m = in our perturbative scalar 
propagators, we have so far ignored the fact that radiative corrections such as the first three 
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diagrams of fig. |T0| will generate a contribution to the mass. And if the mass of a massive 
scalar is treated perturbatively, it will generate infrared divergences, such as depicted by 



fig. [TT], which shows massless propagators connecting a perturbative mass insertion. To be 
at the transition rries = 0, we need to fine-tune the bare mass at this order to satisfy the 



equation of fig. 110. An application of the Feynman rules for these diagrams shows that this 
requires 



'^9^ J fp - j + (higher order) . 



(3.12) 



If we treat mtare perturbatively, it will cancel the radiative contributions to the mass order 
by order in perturbation theory. All we need to do to cancel the infrared divergences in 



our calculation of Aog at this order is to include the additional diagrams of fig. |T^. We will 
henceforth ignore the graphs of figs. ^ and e and the second term of the bare mass (|3.12|) 
above, as these trivially cancel each other. 





(q) (r) 

FIG. 12. Additional NLO diagrams, containing bare mass insertion. 

The graph of fig. ||f gives a contribution to 5\ca of 

-l2Ng'f r / (21-p)'(21 + q)\ 
IZiVS j^J-„j^ ,4(,_p)2(, + q)2 . 

the graphs of figs. @(b,g) and |l^ q give 

(21 + q)2 



(3.13a) 



pq 



l2Ng' fj + (3 - 2Cx)f'f + 3f'f^ 



1^(1 + q)2 



l2Ng^{p^q}; (3.13b) 



and the graphs of figs. §(a,c,h) and |T2]r give 



12Ng' I \3f,f' + (3 - 2Cx)'fj' + ?>flU 



pq 



1 



(3.13c) 
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Now consider the sum of (|3.13| a-c). Setting C\ to 1, we find that the / integral of the sum 
of the integrands converges, as promised.^ 

Having verified that the infrared divergences cancel, it is convenient to proceed by com- 
puting the various terms individually, regulating the I integration with dimensional regular- 
ization. The results are given in Appendix ^. 

The next step is to do the p and q integrations. We do this by first performing the 
integration over the relative angle 9 between p and q. As an example, consider fig. ^ again 
and our result ( p.7|) . By using the angular averages 

' 1 \ 1 , / cos^^ \ p< 



and 



|p + q| 



p> 



|p + q| 



3p> ' 



we obtain 



where 



5>!^ = N(ff :f,,— {i^-ax) 

Jpq p> 



p> = max(p, g) , p< = min(p, q) , x = p</'P> ■ 



(3.14) 



(3.15) 



(3.16) 



Similar results for the rest of the graphs are given in Appendix |B[ The sum of all graphs 
gives 



Jpq Py Jpq 



eff — Ng I ^ pq 

pq P> 



Q + ^x + 3x-^ - (f + ^^'^ + (1 + a;^)"^^^ Sinh'^x . 



6\ 

h{x 

To do one of the integrals easily, write 
h{x) 1 



12 



(3.17) 
(3.18) 



pq 

pq P> 



p dp q dqTpq 

27r* Jo JO p 



1 /•! roo 

- dxx'h{x) dpp\Uflp + fpflp + fphp) 

1 + x)h{x) , 



27r4 JO 
-1 



27r%3iV3c/6 7o X 



1 dx 



l + x , 3 
^l-x) 2 



(3.19) 



^ More accurately, it converges if one ignores logarithmic infrared divergences /^(l/Z^) that vanish by parity. 
The physical regulator — an arbitrarily tiny mass term due to being infinitesimally above the transition — 
respects parity. Dimensional regularization does also. 
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and similarly 



fv fg 



4 



1 dx 



>pq ^ \ An^a^N^g^ Jo x{l - xf 
Combining them gives 

^2 2048 



1 + x^) , 3 So' 



5\ 



cff 



1^9 
27r%3iV2 



where 



K 



1 

X 



; m x H — 

2 



1^1 + x)h{x) 



+ 



6 



{l-xf 



'l + x^ 



l-x) 



■ Inx H X H — x^ 



(3.20) 



(3.21) 



(3.22) 



The contributions of each diagram to k contain logarithmic divergences at small x, but 
the total is integrable. The integral can in principal be done analytically, with the result 
expressed in terms of generalized polylogarithms with arguments like v^, but this seems 
unhelpful enough that we haven't bothered.0 Numerical integration gives 



K = 1.68536... 



(3.23) 



The final results for the tricritical value of X/g"^ is then 



w 2 96 
X g = — 

96 



64k 
37r"2- 



N 



-1 



3.64293iV-^ + 0{N-^) 



(3.24) 



IV. SOME REGULARIZATION ISSUES 

In the last section, we did not bother to introduce a consistent regularization of diver- 
gences diagram by diagram, arguing that that divergences cancel when all diagrams were 



^ The first step is to change variables to y = x + \/l + in the terms involving Sinh ^x, which then 
becomes Iny. This transforms the integral into a product of rational functions and logarithms. 
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summed. This is a potentially dangerous argument and shall later plague us in the SU(2) 
calculation if we do not address it. To see the problem, consider the contributions of the 



form of fig. [T^. The gauge loop integration would have an infrared divergence of the form 

ftm (4.1) 



if the self-energy 11 did not vanish at zero momentum. Fortunately, n(0) = is a consequence 
of Ward identities. The subtlety arises because n(0) does not necessarily vanish diagram 
by diagram. As an example, consider the one-loop contribution to n(0), as shown in fig. |l|. 
Using dimensional regularization for the UV and a small scalar mass for the IR, these 
diagrams give 

If we ignore the UV regularization and instead set (i = 3 and combine the integrands, our 
/ integral is still UV divergent. Our calculation ( [4.1| ) would then end up with ill-defined 
products of IR and UV divergences of the form 

'4- (4.3) 

IR p-^ JUV 



FIG. 13. A class of diagram with IR divergences in the gauge momentum. 
The solution to avoiding this problem is to return to consistently regulating the theory. 



Then rewrite fig. |13 



as 



/ /,'n(o) = / f^[u{p) - n(o)] + / f^u{o) . (4.4) 
jp jp jp 
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The p integration for the first term on the right is now IR convergent, and the calculation of 
n(p) — n(0) is UV convergent diagram by diagram; so we can set d = 3 in this calculation 
and proceed as we did in the previous section. The second term on the right-hand side 
requires full regularization, but we know from Ward identities that n(0) = 0, so we do not 
need to calculate it. The final prescription, then, is to replace Il{p) by Il{p) — n(0) when 



computing diagrams like fig. |13 



The effect of this prescription on the U(l) calculation is that we should add by hand to 
( p. 17] ) a term proportional to 



fpfq 

pq 



(4.5) 



in order to remove any remaining IR/UV divergences in the p/q integrals. However, no such 
term is needed, and so the prescription has no effect on our previous calculation.^ As we 
shall see, however, the prescription will be important to the SU(2) calculation of the next 
section. 



V. SU(2) TRICRITICAL POINT: NEXT-TO-LEADING ORDER 

The SU(2) case is more convenient to treat in Landau gauge than in Feynman gauge. 
This is because there are a host of diagrams, such as fig. 0, which vanish in Landau gauge 
for zero external momentum (because the gluon polarization the external scalars couple to 
is proportional to the gluon four- momentum and does not propagate in Landau gauge). 

We will restrict our attention to the group SU(2) because it has the convenient prop- 
erty that scalar insertions such as shown in fig. |15| are proportional to SabSij. This greatly 
simplifies the analysis of the group factors in diagrams.^ 



^ This isn't an accident. The two loop contribution to ny^^(O) are UV log divergent diagram by diagram. 
For log divergences, cancelations are maintained even when the UV regularization is removed. Consider the 
example of ( |4.2| ) for d—2. 



The quartic interaction |$|'' in (2T) will mean (^J, • $")^, where a is the SU(2) index. As noted in 
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FIG. 14. A diagram which vanishes in Landau gauge. 



a,i 




FIG. 15. Scalar insertion on gauge Hne. 



The Landau gauge results for the Abehan diagrams of fig. ^ are given in Appendix [B 2 . 
For SU(2), they give the contribution 



5A, 



(a-k,q,r) _ ^^^8 
eflf 



•^Pl P> •'fa 



ipq 



27 



(5.1) 



hl{x) — 64 + 256^ 160^^ 128^ ^ + (25^6 + 128^ ^ + 128^^) Siuh X. (5.2) 



In the SU(2) case, 



1 

32 ' 



(5.3) 



and then 



5X. 



(a-k) 
eflf 



214,. .2 



K-19 
TT^ATS ' 



(5.4) 



ref. 1^, the theory with this interaction has a bigger global symmetry than simple flavor IJ{N). The scalar 
sector has an 0(2iV) symmetry. The gauge interactions reduce the symmetry to SU(2)l x Sp(2A^)r, where 
Sp(2A^)r is the A^-flavor generalization of the usual SU(2)r ~ Sp(2)r custodial symmetry of the A^=l case. 
This symmetry is larger than U(A^) and requires that the scalar potential be a function of only the single 
variable • forbidding other possibilities such as • $'^)($jj • €>"). The appearance of such couplings 
severely complicates the generalization of our treatment to general gauge groups. 
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where 



1 dx 

X 



[l + x) 3 

mx + - 

[1-x) 2 

27 



[l + x)hi{x) 



64(1 
-0.052539(1) 



;i + x3), 3 3 2 

(13^ ^"^+2-^+2^ 



(5.5) 



Now consider the non-Abehan graphs, and start with the ghost graph of fig. Doing 
the ghost loop integration, one finds 



(5.6) 



By our prescription for handhng divergences, discussed in the previous section, the second 
term on the right-hand side should be discarded. For SU(2), 



abc cabc 



The figure-eight graph of fig. ||p gives 



6. 



(5.7) 



5X 



(P) _ fabc 4!abc^2 
eS 



jabc jabCgl 



+ / /J/p(!-!cos^^ 



Again, the second term is thrown away by our prescription. 

Finally, we have the graphs of fig. ^1-n). Fig. ^n), for example, is 



(5.8) 



^fabcfabCg.^A I ^^^^ 



32 Jpq 

(2/i - p)^,{2h + g).(2/i - p + q}p I (2/2 - p)-^{2h + g)p(2/2 -p + q)p 

J12 



qih - P)2(li + q) 



/i(l2-p)2(l2 + q)2 



X 



MM 



(5.9) 



where 



r = p + q 



(5.10) 



and 



77 _ Q -fS / / , o f2 f2 . 



pgr 



pJq. 



(5.11) 
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(If one prefers, one can symmetrize the definition of Tpqr with respect to permutations.) We 
have evaluated the scalar integrals by brute force, and various Feynman and momentum inte- 
grals required for this evaluation are tabulated in appendices |A 3| and |A4| . The complicated 
results of the scalar integrations were contracted with the Landau gauge vector propagators 
using a symbolic manipulation program. The result can be most simply expressed when all 
three graphs are combined: 



Jpq 



(5.12) 



where 



^2(p, r) = -{p + q - r){p - q + r){-p + q + r) 



and 



X 



+ 



/i(p + q + r + 2fi) 



{p + q + r)^ Apqr 
{p + q + r + 



32p2g2j.2(^p _|_ g _|_ 



p + q + r + Qp q + Qp r +6g 



2 2 
r 



(5.13) 



(5.14) 



Our prescription for dealing with divergences replaces 



5X^^^^-) ^ fabcfabcgs / Ijr^^^h^ip, g, r) + 3/3g-2 sin^ e„ 

Jpq 

= /'^'7'^''//i-'«2 , 



(5.15) 



where 



i^2 = 7n^i dpdq d{cos 6) p'^q'^J^pqrh2{p, q,r) + 3p^ f '^ sin^ 6 



^71^ Jo 



1 



-0.01573(1) 



(5.16) 



and the value of K2 has been obtained by direct numerical integration.|^ 



The integral can be reduced to a two-dimensional integral by the rewriting p< — xpy and integrating 
analytically over p, giving something too ugly to reproduce here. 
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Putting everything together, the final result for the tricritical value ol \/ is 



A// 



36 
36 



512 



3 " 97r2 
iV"^ + 20.8iV-2 + O(Ar-3) 



4096 , 128 , 2 ^/.r 



37r2 



(5.17) 
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APPENDIX A: USEFUL INTEGRALS 



1. Scalar integrals 



^ij — I Pi Pj 



(Al) 



^(1-Pl)2(l-P2)2 
1 



8Z12 



1 



^(1-Pl)2(l-P2)2(l-P3; 
1 



1 



1 



// P{\ - Pl)2(l - P2)^(l - Pi - V2Y 8P1P2P1 ■ P2 VIPi ~ P2I |Pl 

Using dimensional regularization for the infrared: 

1 



P2I 



1 



//4(l_p)2 
1 











Pi ■ P2 



and rewrite Pi— Ps and P2— Ps in terms of 



qi 



Pl - P3 



q2 



P2 - P3 



\2 ■ 



iPi - Psr IP2 - P3j 

The result is proportional to the bubble integral (|A^), which is trivial to evaluate. 



(A2) 
(A3) 
(A4) 

(A5) 
(A6) 

(A7) 



/^/4(i_p^)2(i_p2)2 8p?pi|pi-p2| 

There is a very simple way to evaluate the triangle integral ([A3|) . First shift 1 — » 1 + ps. 
Then simply change variables again by a conformal inversion: 



(AS) 



(A9) 



The same technique can be used for the integral (|A7|) , which is easiest to implement by 
instead considering the convergent integral 



1 



1 



The inversion relates this to an integral of the form 



(AlO) 
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(1-P1P(1-P2)2 P 



Pi • P2 

8^12 



(All) 



which is easy to evaluate. 



The integral (|A4|) can be reduced to the others by rewriting the numerator in terms of 
denominators as 



1 



1 



/^ + (1-Pl-P2)^-(1-Pl)'-(1-P2)^ 



2pi ■ P2 

It is the only specific case we need of the more general result 

1 



- pi; 



2(l-p2)2(l-p3)2(l-p4)2 



1 



8(^12^34 + ^13^24 + ^14^23) 



1 



1 



1 



+ 

2^12^13^14 ^21^23^24 Z^lZz2ZzA Z41Z42Z43 



(A12) 



(A13) 



2. Angular averages 



1 



1 



1 



|p + q| 

cos 6* 

|p + q| 
1 \ 



cos^ V IP ~ ql IP + ql 



1 

p> 

p< 

2 



Sinh 



-iP< 
P> 



(AM) 
(A15) 
(A16) 



3. Feynman parameter integrals needed for SU(2) scalar loops 



{yzwi + ZXW2 + xywlY 



The Wi are positive scalars. 
2n 



^3/2 [1] 

-^3/2 H 

•^3/2 [a;^] 



W1W2W3 



2tt 



W2W3{Wi + W2 + W3) 



TT 



W2W3 



Wi 



{Wi + W2+ W3) {Wi + W2+ 



(A17) 



(A18) 
(A19) 

(A20) 
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^3/2 M = 



TT 



•^3/2 [■ 



X 



•^1/2[1] 
^1/2 W 



TT 



TT 



+ 



3wi 



+ 



2wl 



4wa 



{W1+W2 + W3) {Wi+W2 + W3y (Wi+W2 + W3)3 
1 2Wi 



+ 



{Wi +W2+ Wa)^ (wi + W2 + Ws)^ 



TT 



2(^1 +W2 + Ws)^ 
71 

W1 + W2 + W3 
71 

4 



+ 



Wi 



(A21) 

(A22) 

(A23) 
(A24) 
(A25) 

(A26) 



{W1+W2 + W3) (wi + W2 + Ws)'^ _ 

The ^3/2 results can be obtained from the J^i/2 results by differentiating with respect to the 
Wi and using x + y + z — 1. 



4. Scalar integrals needed for SU(2) 



J = ^i2P3 + 2;23Pl + 2;3lP2 

K = pi + p2 + Pa 



(A27) 
(A28) 



Ji 



I (I - pi)2(l - P2)^(l - P3)^ 
(I-Pl)2(l-P2)2(l-P3)^ 



8^12^23^31(^12 + ^23 + ^31) 

5. 

+ 



Jij 



JiJj 



(A29) 



16(2;i2 + 2:23 + 2:31) 162;i22;23^3l(^12 + Z23 + ^3l)^ 
Zl2P3iP3j + Z23PliPlj + Z3iP2iP2j 



liljlk 

'.(l-pi)2(l-p2)2(l-p3)2 

J I Jj J }^ 



^ l'0Zi2Z23Z3i{Zi2 + Z23 + Z^l) 

JiSjk + JjSki + Jk^ij Ki5jk + Kj6ki + Kk6ij 



(A30) 



64(2;i2 + ^23 + ^31)^ 64(2;i2 + ^23 + ^31) 

{Qzf2P3iP3jP3k + lSzi2Zi3P2iP2jP3k) 



32zi2Z23Z3i{zi2 + 2^23 + 2:31)^ (i4:Zi2Z23Z^i{zi2 + 2:23 + 2:31)^ 

{Qzi2P3iP3jP3k) 



(A31) 



()4:Zi2Z23Z^l{zi2 + 2:23 + 2:31) 

The angle brackets in the last expression denote averaging over all permutations of k) 
and all permutations of (pi,P2,P2)- 
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APPENDIX B: SUMMARY OF DIAGRAMS 



1. U(l) case: Feynman gauge 



6X 



(a,c,h,r) 
eff 



6\ 



(f) 

eff 



6\ 



(i) 
eff 



6X 



(k) 
eff 



eq. (3.13c) 



pq 



3fJ^ + (3 - 2C,yf^f^' + 3f^f, 



pJq 



IpJq 



TT^A 



24 



TT^A 



^A^S.'^'"^) = eq. { ^Jm 



pq 



Uf^ + (3 - 2C,)/pV; + 



'6 ^ 12 



q TT^A 



pq 



pq 



1 12 

— (3 + 3x-i) + — - 



eq. (|3.13a|) 

Ng"" I Tpq 

Jpq 



6 6 



6 



3cos^ ^3(j92 + g2) 



6 



V 1 |p + qI 2|p + q| P9|p + q| TT^A 



v> 



Ng' / 

Jpq 

eq. (P) 
= iV// J-p,— (-9) 

Jpq p> 

a2 = eqs. e 



^ ' 12 - 3a;-i + f a;) - 



iV/ / J'pq— {IS -Ax) 

Jpq p> 



(21 - p) ■ (21 - p + 2q) (21 + q) • (21 - 2p + q) 
/2(l-p)2(l + q)2(l_p + q)2 



Jpq 
pq Jl 



6 6 9(p2 + g2) / 
p q 2pq \IP + ^ 



1 1 

+ 



IP - q| 



+ 



3(2p^ + 5p2g2 + 2q^) 3p ■ q 



1 



(Bl) 



(B2) 



(B3) 



(B4) 



(B5) 



1 



Ng , . pq 



pq 'P> 



Apqp ■ q 

-6 + 3x~^ + 7x 



pq 



|p + q| |p-q| 



(B6) 



The graphs of figs. ^ and e and the second term in the mass counter-term ^.12| have been 
ignored above, as they trivially cancel each other. 
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The infrared divergences of ( pi|) through (|B3|) were regulated with dimensional regular- 
ization, but, for the sake of making cancelations explicit, we have put the linear divergences 
back in by hand by writing 

where A is an infrared momentum cut-off. 



2. Abelian graphs: Landau gauge 



Six) = 



1, U(l) theory 
X, SU(2) theory 



fB8) 



r\(a-d,g,h,q,r) 



6X 



6X 



6\ 



(f) 
cflf 

(i) 
eflf 

(j) 
eflf 



(k) 
cff 



Jpq q 

{is)N9' f -F,,-^ (-8 + 4x - 

^ ^ Jpq P> ^ 

Ng' f ^p,— (-4- 

Jva ^ 

128 / ^9 

' (-lis) ^9 



^ ( 128 



Jpq p> 
1 

PI 

pq P> 



4 2 
5 



2 2 
5-^ 



4x + |a;^ 



T 

■Jpq 

pq P> 



^x 



f + 3x-^ + 3x2) ^ ^2^-1/2 ginh-ix 



(B9) 
(BIO) 
(Bll) 
(B12) 

(B13) 



For non- Abelian graphs, see eqs. ( ^.6|) , (|5.8|) and ( ^.12|) . 
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